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ABSTRACT 

O- and B-type stars are often found in binary systems, but the low binary mass-ratio regime is 
relatively unexplored due to observational difficulties. Binary systems with low mass-ratios may have 
formed through fragmentation of the circumstellar disk rather than molecular cloud core fragmen- 
tation. We describe a new technique sensitive to G- and K-type companions to early B stars, a 
mass-ratio of roughly 0.1, using high- resolution, high signal-to- noise spectra. We apply this technique 
to a sample of archived VLT/CRIRES observations of nearby B-stars in the CO bandhead near 2300 
nm. While there are no unambiguous binary detections in our sample, we identify HIP 92855 and 
HIP 26713 as binary candidates warranting follow-up observations. We use our non-detections to 
determine upper limits to the frequency of FGK stars orbiting early B-type primaries. 
Subject headings: binaries: spectroscopic — stars: early-type — stars: formation 



1. INTRODUCTION 

O- and B-ty pe stars are often fo und in binary or mul- 
tiple systems: IMason et al.l ()2009l ) estimate that at least 
57% of 0-stars are in spectroscopic multiple systems, 
and at least 75% are in any type of binary or multiple 
system. Yet the multiplicity fraction of high-mass stars 
may be underestimated du e to the difficulty of d etecting 
low- mass secondary stars (jSana fc Evand |20TTI ) . While 
the mass-ratio distribution is reasonably well known for 
high-mass binaries with mass-ratio q = M^/Mp > 0.2, 
there are almost no constraints for low mass-ratio bina- 
ries. However, binaries of low mass-ratio are important 
probes of star formation since they may have formed in 
a different way than approximately equal-mass binaries. 
Here we define low mass-ratios to be those with q < 0.2, 
where Ms is the mass of the secondary (lower mass binary 
component), and Alp is the mass of the primary (higher 
mass binary component). We also use the term "low- 
mass" to describe star with M < IMq, and "high- mass" 
to describe stars with M > 5Mq 

1.1. Binary Formation in High- Mass Stars 

With such a high fraction of high-mass stars found 
in binary or multiple systems, any theory of high-mass 
star formation should be able to explain the high bi- 
nary formation rate. There are several mechanisms by 
which binary s tars m ay form: fission ()Lvttletonl 119531 : 
iLebovitj 11974 I1984D . in which a molecular core begins 
spinning fast enough as it collapses that it splits into two 
stars ; core fragmentation (see e.g. iBoss &: Bodenheimeri 
119791 : iBossI 119861 iBate et al.l[l995f ). in which a collaps- 
ing core develops two or more overdensities which then 
begin collapsing separately; and disk fragmentation (see 
e.g . Krattcr & Matzner 2006; Stamatcllos & Whitworthi 
[201 ll) ■ in which the circumstellar disk surrounding the 
primary star becomes gravitationally unstable and cre- 
ates a secondary star. While the fission scenario was 
once thought to be important, it has since fallen out of 
favor because the viscous dissipation timescale, which 
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would drive a spinning body towards fission, is much 
longer than the core collapse timescale ( Tohline ,200^ 
and because hydrodynamic simulations fail to cause the 
rotating core to actually split rather than just deform 
(jTohline fc Durisenl[200ll ). Both core and disk fragmen- 
tation are still thought to be viable binary formation 
mechanisms. It is likely that both mechanisms play a 
role in shaping the binary mass-ratio and separation dis- 
tributions. In the formation of higher-order multiples, 
it is even possible that both mechanisms operate in the 
same system. 

The primary method of forming binary systems is 
thought to be core fragmentation. As a molecular cloud 
begins isothermally collapsing, its density increases, 
causing the Jeans mass to decrease. Thus, an initially 
Jeans-mass collapsing core can fragment into smaller ob- 
jects. Core fragmentation will initially yield binaries with 
separations lOAU < a < lOOOAU, which may move 
closer by interacting with the surrounding gas, a cir- 
cumbinary disk, or through dynamic al interactions with 
other nearby stars (|Bate et al.l 120021 ) . For ~ IM© pri- 
mary stars, the observed mass-ratio distribution is well 
fit by assuming both components of the binary are cho- 
sen randomly from the stellar initial mass function, and 
later evo lve through accr etion and dynamical interac- 
tions ( Kroupa et al.|[2003f ). Assuming independent com- 
ponent masses chosen from the IMF, a binary with a 
lOAf0 primary would most often have a 0.1 Mq sec- 
ondary, giving an initial mass ratio near q = 0.01. The 
unmodified mass-ratio distribution of high-mass binaries 
would therefore strongly favor low mass ratios. How- 
ever, accretion of high specific angular-momentum gas 
from either the collapsing molecular core or the circumbi- 
nary disk will preferentially be captured by the lower- 
mass companion, driving the binary mass-r atio toward 
unity and decreasing the orbital separation (|Batel[2000l : 
iBonnell fc Batd 120051 ). Additionally, dynamical interac- 
tions tend to replace low-mass binary companions with 
higher-mass ones, or to kick the lower-mass component 
out to a wide orbit and create a hierarchical triple sys- 
tem. Over time, these processes tend to create high-mass 
binary systems with nearly equal masses and small sep- 
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arations (|Bate et al.ll2002t) . Dynamical interactions are 
most important in dense environments where the prob- 
abihty of stellar encounters is high. Therefore, they are 
probably more important in dense OB star clusters than 
in the much looser OB associations, where many B stars 
are found. 

Another potentially important way to form binary 
systems is disk instability (see e.g. iKratter fc Matzneii 
[20061 IStamatellos fc Whitworthll2Cillir In this scenario, 
the fragment forms in an unstable circumstellar disk 
with an initial separation of ~ lOOAU and initial mass- 
ratios nea r q ^ 0.03, similar to t he core fragmenta- 
tion case (|Kratter fc Matzneii 120061 ). The final mass- 
ratio is expected to rise, but not as significantly as in 
the core fragmentation scenario in which the fragment 
can form sooner. Typical mass-ratios near q ^ 0.1 are 
expected, though the picture is far from complete. A 
se mi-analytical tr e atmen t of embedded protostellar disks 
by IKratter et al.l ()2008fl finds that massive stars with 
M > 2Mq maintain 0.01 — 0.1 in o rbitin g fragments 
after about 2 Myr. iKrumholz et al.l ()2007f) simulate a 
IOOM0 collapsing core for a much shorter time (20 kyr), 
but also find that the disk fragments a nd that the final 
fragment mass ratio is g ~ 0.1. However . IKrumholz et al.l 
(|2009i) simulate the same mass core but start it with 
a slow solid body rotation instead of a turbulent ve- 
locity field and run the model for about twice as long 
(57 kyr), and find that it leads to a very massive bi- 
nary (M l -I- M2 =_70Mq) with a mass-ratio q = 0.7. 
Work by IClarkd (|2009[) indicates that as mass is trans- 
ported inwards onto the star, the outer disk can be- 
come unstable at late times. This instability can lead 
to a delayed disk fragmentation, with a fragment mass- 
ratio in the r ange 0.1 < q < 0.5. The simulations by 
IClarkd (|20?)9h were done for a ~ IMq primary star, 
and the delayed fragmentation occured after about 10^ 
years. Delayed fragmentation may not be possible in 
disks surrounding high-mass stars, as the time at which 
fragmentation occurs is comparable to the disk dispersal 
timescale (Klahr & Brandncr 2006). However, if it does 
occur, the similar fragmentation and dispersal timescales 
suggest that the fragment would not undergo significant 
accretion or migration and would leave a wide binary 
(a w lOOAU) with mass-ratio in the range 0.1 < q < 0.5. 

Unfortunately, there are no true binary population 
synthethis simulations for high-mass binary systems 
formed by either mechanism discussed above. The lack of 
population synthethis models is driven by computational 
issues; a collapsing cloud that reproduces the stellar mass 
function and generates enough high-mass stars to mean- 
ingfully analyze the binary statistics would have to be 
very massive and therefore difficult to simulate. Disk 
fragmentation simulations often either stop the simula- 
tion once a fragmen t forms rather than follow its mass ac- 
cretion history (e.g. lBossll201 It IKrumholz et al.ll2007l ). or 
lack high enough resolution to follow the secondary very 
near the star (e.g. BonncU & Bate 2005). Fortunately, 
such models may be in the near future. Realistic simula- 
tions of massive collapsing molecular clouds have begun 
appearing that can meaning fully discus s the multiplicity 
of low -mass binary systems (jBatd [20121 : IKrumholz et aJ] 
I2OI2I ). These simulations reproduce the observed in- 
crease in multiplicity fraction with primary star mass. 



but do not yet generate enough high-mass binary sys- 
tem to compare the parameter distributions to observa- 
tions. Despite the current lack of models, we can draw 
the general conclusion that disk fragmentation tends to 
produce lower-mass companions than core fragmenta- 
tion. For this reason, probing the low mass-ratio regime 
can provide information on the relative importance of 
both scenarios in forming high-mass binary systems, and 
may help constrain models once computational power in- 
creases. 

In addition to binary star formation, disk instabil- 
ity is often invoked as a way to form planets of a few 
Jupiter masses orbiting ~ IMq stars. While the mas- 
sive star formation process as a whole may not sim- 
ply be a scaled up versio n of low- mass star formation 
ijZinnecker fc Yorkd[200l . the process of disk fragmen- 
tation may be. One expects that disks around high- 
mass stars, with correspondingly higher accretion rates 
and more mass . frag:ment more often than disks around 
low- mass stars (lBosg|201 11120061 : iDodson- Robinson et al.l 
i200a IKratter fc Matzneiil2006[) . Thus, if disk fragmenta- 
tion plays an important role in high-mass star formation, 
it may also play a role in low-mass star formation by cre- 
ating ^ lOMjup planets and substellar companions. 

1.2. Observing low mass-ratio binaries 

Detection of OB-star binaries with mass-ratio q ^ 0.1 
or lower is very difficult, since the ratio of the secondary 
flux Fs to the primary flux Fp is Fs/Fp ~ 10"'^ or lower 
in the V-band. Imaging surveys can detect such contrast 
ratios for wide orbits, but lose se nsitivity as the separa - 
tion decreases below about 1" (e.g. lMafz Apellanizll2010f) . 
Spectroscopic binary surveys do well for short-period sys- 
tems where a full orbit can be mapped in a reasonable 
amount of time, but lose sensitivity for periods greater 
than about one year (e.g. iSana et al] 120091 : lEvans et al.l 
l2010f ). However, low-mass companions (g < 0.2), which 
induce a small reflex motion on the primary, are very 
difficult to find with traditional spectroscopic surveys. 

One method to find low-mass companions to late B- 
type primaries is to search for high x-ray emission. 
Stars later than about B3 are no t expected to hav e 
strong enough winds to emit X-rays (jGagne et al.lll997D . 
and stars of earlier type than about A7 do not have 
a radiative-convective boundary that can drive a mag- 
netic d ynamo and create an X-ray generating corona 
(jSchm itt 1997). Stars in between spectral types B4 and 
A7 with strong X-ray emission are thought to have young 
low-mass companions, because the luminosity and X- 
ray spectral e nergy distribution is similar to observed 
T-Tauri stars ([Huilamo et al.|[200l . lEvans et al.l poTTI ) 
use this fact to search fo low-mass companions to late 
B-stars in the open cluster Trumpler 16. They find a 
significant number of companions, and set the multiplic- 
ity fraction at 39%. This value is a lower limit, but the 
authors believe that the true value is not much above 
39%. Unfortunately, X-ray imaging is not effective for 
primary star spectral types earl ier than B3, which are 
also strong x-ray emitters (Gagne et al.l Il997l ) and will 
drown out any companions. 

In this paper, we introduce a technique that is sensi- 
tive to young binary systems with secondary tempera- 
tures 4000 K < Toff < 6000K. For early B-type primaries 
with ages ~ 15 Myr, these temperatures correspond to 



Detection of Low Mass-Ratio Stellar Binary Systems 



3 



Detection Method 



^ 0.985 
0.980 





2324 2326 2328 




Velocity (km/s) 



Fig. 1. — This figure illustrates the approximate fiux ratio limit to the detection method outlined in section|2] Top panel: Residuals after 
telluric correction (see section |4]l for chip 2 of HIP 80582 are in black, with an atmosphere model for an 0.9Mq star at 50.1 Myr below it 
in red. The flux ratio at this age is Fs/Fp = 0.0092. Middle panel: The scaled model spectrum was added to the telluric residuals, and 
then the sum was cross-correlated with the model. Despite the signal being significantly below the noise level, the star was detected at a 
high significance. The y axis, in units of the standard deviation of the cross-correlation function, shows that the significance of the peak is 
over 4cr. Bottom panel: Same as the middle panel, but the model spectrum was added to the residuals with a 50 km s~^ velocity offset. 



mass-ratios q w 0.05 — 0.3, right where we expe ct t o 
see binaries formed by disk instability (see section [TTT|) . 
Rather than attempting to detect the reflex motion of 
the parent star as in exoplanet searches and SBl bina- 
ries, we attempt to directly detect the spectrum of the 
young low-mass companion using high signal-to-noise, 
high-resolution data. There is a multitude of archived 
B-star observations in the near-infrared, where they are 
used as telluric standard stars to remove the absorption 
spectrum of the Earth's atmosphere (telluric lines) . This 
method is equally sensitive to all separations within the 
point spread function, which is dominated by the seeing 
since the adaptive optics are not used in telluric standard 
star observations. A typical seeing of ~ 0.8" corresponds 
to up to ~ 900 AU for targets within a few kpc. In this 
paper, we describe a search for young F5-K9 type com- 
panions in archived VLT/CRIRES spectra of 34 early 
B-type stars. 

In section [2] we describe our detection method in more 
detail. Section [3] describes the B-star sample we use in 
this work. Section 2] contains the data reduction and 
telluric correction methods. We summarize our results 
in section [5l We examine the completeness of our sample 
in section [5] and put limits on the multiplicity fraction as 
a function of mass-ratio in section [T) Finally, we present 
our conclusions about the prevalence of low mass-ratio 
companions to early B-type stars and discuss how our 
results constrain star formation mechanisms in Section 

E 

2. DIRECT SPECTRAL DETECTION METHOD 

We describe here our method to detect the emission 
from an approximately solar-mass star orbiting an early 
B-type star, which we will hereafter call the direct spec- 
tral detection method. The basis of this method is to 
cross-correlate a high signal-to-noise ratio B-star spec- 
trum with a synthetic F,G, or K star spectrum. If a 



low-mass star with such a spectrum is orbiting the B- 
star, we expect to find a peak in the cross-correlation 
function at the radial velocity corresponding to the low- 
mass star's motion. A peak in the cross-correlation func- 
tion should appear even if the flux from the low-mass 
star is comparable to or even slightly less than the noise 
level in the spectrum. Figure [T] illustrates the approx- 
imately limiting case for the flux ratio. The top panel 
shows a fully reduced CRIRES spectrum of HIP 108975 
(see sectional) with a model spectrum for an O.9M0 star 
at a realistic flux ratio below it. We used evolution- 
ary tracks published by iLandin et al.l ()2008() to evolve 
the secondary star to 50.1 Myr, the age of the system 
(jTetzlaff et al.|[20Tot) . in order to determine the flux ra- 
tio between the primary and secondary. The secondary 
star model was then added to the telluric-corrected B- 
star spectrum at two different velocities. It is clear that 
the model spectrum has an amplitude much smaller than 
the noise. Nonetheless, the bottom two panels show that 
a cross-correlation will have a peak with high significance 
at the velocity of the secondary star. 

A careful choice of the wavelength region is critical for 
the direct spectral detection method. First, we want a 
wavelength region where the B-star spectrum is mostly 
continuum (i.e. very few spectral lines). Since B-stars 
have few spectral lines, it is easy to find a such a spectral 
region. Secondly, we want a region where the low-mass 
star would have many closely spaced, strong lines. The 
more lines there are in the low-mass star, the stronger the 
peak will be in the cross-correlation function. Finally, we 
want a spectral region where the flux ratio between the 
low-mass and the high-mass star is maximized. It is not 
helpful to go much redder than a few microns for com- 
panions with T > 4000K, because both the high-mass 
and the low-mass star arc firmly in the Rayleigh- Jeans 
limit by this point, where the flux ratio is approximately 
constant. For this project, we choose wavelengths from 
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Fig. 2. — Top panel: The telluric spectrum (absorption due to Earth's atmosphere) in the wavelength range from 2290 — 2400 nm. Most 
of the lines are from CH4, with a few H2O lines appearing in the right half. Bottom panel: The model spectrum of a 5500 K star with 
log((;) = 4.0 and solar metallicity. Note that the line density of telluric lines is comparable to or greater than that of the star model, and 
many of the telluric lines are stronger than the stellar lines. 



2300 - 2400 nm, which is the CO Ai/ = - 2 bandhead 
in the low-mass star. 

There is both a lower and upper mass detection limit. 
Secondary stars that are too cool will be too faint, and 
any signal will be lost in the noise. Additionally, a more 
massive (and hotter) primary star will decrease the flux 
ratio, and push the lower mass limit up. On the other 
end, secondary stars that are too hot will dissociate CO, 
destroying the bandhead that we are looking for. The 
temperature and size of the secondary star will depend 
on its age as well as its mass since it will still be evolving 
towards the Main Sequence during the lifetime of the B- 
star. The exact mass sensitivity will thus depend on the 
age, primary star mass, and signal-to-noise ratio of the 
system being observed. 

The detector resolution is also important for the di- 
rect spectral detection method. Deeper lines, providing 
more contrast from the continuum, are easier to detect 
than broad, shallow lines. In addition, narrow spectral 
lines will result in a stronger, narrower peak in the cross- 
correlation function, which is most sensitive to the steep 
line edges. Therefore, we want the spectral lines in the 
low-mass companion to be as deep and narrow as possi- 
ble. The intrinsic width of CO bandhead lines is roughly 
5-7 km s~^. In order for the observed line width to be 
this small, we need the resolution of the instrument to 
be i? = A/AA > 50000. 

There are two main difficulties with the direct spectral 
detection method: telluric line removal and the low flux 
ratio between the primary and secondary star. Figure 
|2| shows the transmittance through the Earth's atmo- 
sphere (the telluric spectrum) in the wavelength range 
we are interested in. Most of the spectral lines are from 
methane, with a few deep water lines towards the red end 
of the range shown (See section |4| for details on the tel- 



luric line removal). The low flux ratio makes the telluric 
contamination especially troublesome, since the telluric 
lines are stronger than the lines in the companion star 
spectrum. The flux ratio of Fg/Fp ~ 10"^ effectively 
sets a lower limit on the signal-to-noise ratio for which 
the direct spectral detection method is possible. Any flux 
coming from a low-mass star will be completely buried in 
the Poisson noise for spectra with SNR ^ 100. Removal 
of the telluric contamination will add more noise to the 
spectrum, so a spectrum should have SNR of a few hun- 
dred before telluric line removal to have a good chance 
of detecting a companion. 

3. STAR SAMPLE 

B-type stars are commonly used in the near-IR as tel- 
luric standard stars. Astronomers will observe their sci- 
ence targets, and then move to a B-type star. Since 
B-type stars have few spectral lines relative to cooler 
stars, most of the observed spectral lines will be from the 
absorption of Earth's atmosphere (telluric absorption). 
Therefore, these stars provide an empirical estimate of 
the telluric spectrum; division of the science spectrum 
by the normalized standard star spectrum will mostly 
remove the telluric lines. 

Since B-type stars are commonly used as above, there 
are many high signal-to-noise ratio (S/N > 100) obser- 
vations of such stars in archived data. We used the 
VLT/CRIRES archive in this project. CRIRES is a high 
resolution {R = A/AA « 100000) infrared spectrograph 
on the VLT at Paranal Observatory. The detector con- 
sists of four 1024x512 ccd chips that are mosaiced end- 
to-end, and the spectrum falls across them. There are 
several wavelength settings available, which determine 
what parts of the spectrum fall on each chip. For wave- 
length settings in the CO bandhead near 2300 nm, each 



Detection of Low Mass-Ratio Stellar Binary Systems 



5 



chip will hold rougly 10 nm of spectrum with roughly 1-2 
nm gaps between the chips. 

To generate the sample, we started with all single, 
main sequence B0-B5 stars with CRIRES observations 
from 2300-2400 nm. We then excluded any s hell stars, 
which have circumstellar disks (jPorter fc Rivinius 2003) 
that may create false positives. Table [T] shows the com- 
plete sample used in this project. The spectral types 
and a ges were obtain ed from a catalog of nearby young 
stars (|Tetzlafr etaD l2010|)^ The one exception is HIP 
97611, which was not in [Stzlaff ct al. (2010). For this 
star, the age was taken from fWestinI (Il985i) and the spec- 
tral type from the Simbad database! ■ The distances to 
all stars were determined from parallaxes given in the 
Simbad database. The maximum separation column es- 
timates the approximate maximum separation of the bi- 
nary orbit we are sensitive to, assuming a seeing of 0.8" 
which is typical of Paranal Observatory. The median of 
the maximum separations to which we are sensitive is 
124 AU. The final column gives the number of distinct 
observations of the star. We count all nodding positions 
taken on a given night with the same detector wavelength 
setting as one observation. 

4. DATA REDUCTION AND TELLURIC CORRECTION 

The data reduction was done using standard meth- 
ods in IRAlfl. All observations were taken in an AB or 
ABBA nodding pattern. For each set of AB nods, A-B 
and B-A frames were made to remove any atmospheric 
emission lines and dark current. The resulting difference 
images were then treated to a quadratic nonlinearity cor- 
rection, using coefficients made available by the CRIRES 
teamThe corrected frames were then divided by a nor- 
malized flat-field. Due to the slit curvature, the spectrum 
can shift by up to a pixel in the dispersion direction be- 
tween the A and B nod positions. Therefore, combining 
the 2D frames before extraction can reduce the spectral 
resolution and affect the line shapes. For this reason, 
we combined the nodding positions only after the wave- 
length calibration and telluric correction. Each nod po- 
sition was extracted using the optimal algorithm in the 
apall task in IRAF. The spectra were wavelength cali- 
brated using a model telluric spectrum g enerated with 
the a tmospheric modeling code LBLRTM IjClough et al.l 
[20051) . 

For telluric correction, w e used a similar procedure to 
the one outlined by Seifah rFet al.l (2010). The atmo- 
sphere modeling code LBLRTM was used to generate a 
synthetic telluric absorption spectrum. The abundances 
of water, methane, and carbon monoxide were fit using a 
Python implementation of a Levenberg-Marquardt non- 
linear least squares fitting algorithm. The Levenberg- 
Marquardt fit also refined the wavelength solution to the 
telluric model, fit the continuum, and fit the resolution of 
the spectrograph with a Gaussian profile. The FWHM of 
the profile was the only free parameter in the resolution 
fit. 

The LBLRTM code expects a model atmosphere, 
which contains the temperature, pressure, and abun- 

^ |http : // simba d . u- strasbg ■ f r/simbad/ 1 

^ IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



dance of 30 molecules as a function of atmospheric 
height. For the majority of molecular species, we used a 
mid-latitude nighttime MIPA^ profile, which provides 
the temperature, pressure, and abundances of various 
molecules in 1 km intervals from sea level to 120 km. 
The low-altitude ( z < 30 km) te mperature, pressure 
(jKerber et al.ll2010[ ). and humidity (jChacon et al.ll2010D 
profiles were obtained from radiosonde data taken from 
Paranal Observatory. 

The LBLRTM atmospheric modeling code comes with 
a molecular line list based on the HITRAN 2008 database 
(jRothman et al.ll2009t ) , with a few molecules individually 
updated. Since none of these updates were relevant for 
the wavelength range from 2300 — 2400 nm, we in essence 
used the stock HITRAN 2008 database. However, in the 
process of modeling, we found several water and methane 
lines that were consistently under- or over-fit. For these 
cases, we manually adjusted the line strengths in the 
database. The line strengths were fit visually and should 
not be considered rigourous new line strengths. Table [2] 
summarizes these changes. 

In some of the 2007 data, the first chip was not well 
illuminated by the flatfield lamp. This introduced an 
unphysical continuum shape in the data and made the 
resulting model fit very poor. For these cases, we ig- 
nored the first chip in further analysis. In addition, the 
fourth chip has several bad pixels on the left edge and a 
streak down the middle. None of the telluric model fits 
were very good on this chip, and so we have ignored it 
completely in our analysis. 

After the observed spectrum was fit, we found that the 
residuals still contained large spikes, even on the good 
detector chips. These spikes can come from a variety 
of sources. For the deepest lines, simple Poisson noise 
can create large residuals when dividing by the telluric 
model. In addition, a poorly fit continuum may cause 
the model to over- or under-estimate the abundance of a 
given molecule. This can be especially troublesome for 
water lines, for which only a few exist in the wavelength 
region we are investigating. If a strong water line is near 
the edge of the chip, where the continuum is usually least 
certain, the best-fit water abundance may be skewed and 
cause none of the water lines to be well fit. We do know 
that large residuals are not coming from the spectrum of 
a low-mass star, due to the expected flux ratio between 
the primary and secondary, Fg/Fp ^ 10^^. Any residuals 
with amplitude greater than 1% of the continuum level 
come from uncorrected telluric lines, cosmic rays, or bad 
pixels. 

In order to minimize these spikes, we performed a sec- 
ond fit to any residuals significantly above the contin- 
uum noise level. To make sure we were not fitting away 
any low-mass star lines, we only corrected spikes whose 
amplitude was greater than 5% of the continuum level. 
In this second fit, we first attempted to fit a Gaussian 
to each spike. If the spike was well fit by a Gaussian, 
we divided the residuals by the fit. If not, we simply 
masked out the line core, so that it would not affect the 
cross-correlation in later analysis (see section [2]). Fig- 
ure [3] shows the steps involved in the telluric correction. 
Notice that the secondary correction removes the large 
residuals, while leaving the rest of the spectrum unaf- 

^ |http : //www-atm . physics . ox ■ ac. uk/RFM/atm/ 1 
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Fig. 3. — The telluric correction steps for chip three of CRIRES wavelength setting X^ef = 2329.3 nm. Top panel: Normalized spectrum 
(black), with the best-fit telluric model (red). Middle panel: Residuals after dividing the observed spectrum by the telluric model. Note 
the large spikes near 2328.5, 2331, and 2335 nm. Bottom panel: Correction after fitting the large residuals to Gaussians. 

line removal in a multi-stage process. 

The first rejection stage was done by identifying peaks 
in the cross-correlation caused by telluric residuals. To 
do this, we cross-correlated a spectrum uncorrected for 
telluric absorption with the same suite of model atmo- 
spheres as we used for the corrected spectra (see above). 
We did these cross-correlations for one observation of 
each wavelength setting. The cross-correlation of un- 
corrected spectra with model-atmosphere spectra gen- 
erated a series of cross-correlations with peaks arising 
exclusively from telluric lines. We visually compared all 
of the binary candidate signals with these telluric cross- 
correlations. If the dominant cross-correlation peak was 
at the same velocity and had a similar width as a peak 
in the telluric cross-correlation function corresponding 
to the same wavelength setting and secondary model 
temperature, we assumed that the peak was caused by 
incomplete telluric removal and rejected the candidate. 
There were several cross-correlations with peaks at the 
same location as a telluric peak, but with a different 
width. In these cases, we marked the candidate as prob- 
ably coming from incomplete telluric correction, but did 
not reject the candidate. 

Next, we determined whether the signal-to- noise ratio 
and telluric line removal in a given observation would al- 
low us to detect the candidate companion star. To do 
this, we added a model atmosphere with the same tem- 
perature as the candidate to the telluric-corrected spec- 
trum at 17 different radial velocities ranging from -400 
to 400 km s~^. We do not expect to see any peaks from 
real companions with \v\ > 400 km s^^, the approximate 
radial velocity of a IMq star orbiting a IGA/0 star such 
that the stellar surfaces are in contact. The flux ratio 
of the model atmosphere to the primary was obtained 
by in terpolating pre- main -sequence evolutionary tracks 
from iLandin et al.l (|2008[ ) at the age of the system, as 
well at ageicTage- The primary star ages for our sample 
are given in Table [TJ We cross-correlated each of these 
semi-synthetic spectra against the model spectrum; if the 
largest peak was at the correct velocity, we counted the 
star as detected. If the star was not detected in the 
sensitivity analysis at least 50% of the time, we rejected 
the candidate. The significance of the correct peak in the 



fected. 

The telluric correction described above usually reduced 
any telluric lines to near the Poisson noise level in the 
spectrum, which is the best a fitting routine can do. To 
search for any systematic errors in the telluric correction, 
we added all spectra of each wavelength setting together 
to make a series of master telluric residual spectra. These 
master spectra had less random noise than any individual 
observation, and therefore we were immediately able to 
see whether some telluric lines are systematically under- 
or over-fit. We found that for wavelength settings with 
at least ten spectra in our sample, dividing the corrected 
spectra by this systematic telluric residual template in- 
creased the sensitivity to companions. We did not make 
this final correction for wavelength settings with fewer 
than ten individual spectra in our sample. 

5. RESULTS 

Each telluric residual spectrum was cross-correlated 
against a suite of model atmospheres generat ed by 
the P hoenix stellar atmosphere code (jHauschildt et al.l 
IIQQQT ). All model spectra had solar metallicity. The ef- 
fective temperatures ranged from 3000 — 7200 K, in lOOK 
intervals. We used several surface gravities based on 
the stellar temperature. For the model secondary stars 
with 3000 < Toff < 3600, which would have to be very 
young (and large) to be detectable, we used a log(5) = 
3.5. For the secondaries with 3600 < Toff < 6500, we 
used \og{g) = 4.0. Finally, we used \og{g) = 4.5 for 
Tcff > 6500, which can be detected closer to the Main 
Sequence. We found that the surface gravity has only a 
very small effect on the cross-correlation, which is more 
sensitive to the line position than its precise width or 
depth. We compiled a list of all cross-correlations that 
show a single peak with at least 3cr significance. For a 
given telluric-corrected residual spectrum, several model 
atmospheres may generate a significant peak at the same 
velocity. This is because the model spectra of two stars 
differing by only a few hundred kelvin are not very differ- 
ent. To keep from counting peaks twice, we only counted 
the cross-correlation that resulted in the most significant 
peak at a given velocity. We then attempted to reject 
spurious peaks caused by the noise or incomplete telluric 
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Fig. 4. — Cross-correlation for HIP 26713, using a 5600 K star model spectrum as template. The y-axis is in units of the standard 
deviation of the cross-correlation function. The peak is very near the maximum velocity of |wmax| = 256 km s~^, assuming a circular 
orbit (see equation [TJ . The likelihood of observing the system nearly edge on and at a quadrature point, so that li^l ~ |i'max|, is p 0.01. 
However, it is possible that the system has an eccentric orbit, effectively increasing |»;max|. 

able to reject. For both of these, we checked what other 
observations the candidate star had within the CO band- 
head spectral region (2300 — 2400 nm). The analysis of 
each of these stars is done separately below. 



cross-correlation function can vary greatly, depending on 
where the stellar spectrum falls in relation to the telluric 
line residuals. Therefore, we cannot usually reject a peak 
based solely on its significance. 

We then visually inspected the remaining candidate 
cross-correlations, picking out those with a single domi- 
nant peak with [iirl < Wmax where Wmax is the maximum 
possible radial velocity for a star of temperature Tgec to 
be orbiting a hotter star of temperature Tprim with a 
semimajor axis a. Assuming a circular orbit, i^max is 
given by 



'2G(MpHm + M,ec) T,, 



R, 



■prim 



T ■ 

pri 



(1) 



prnn 



where Afprim and Msec are the masses of the primary 
and secondary stars, respectively, and i?prim is the ra- 
di us of the primary sta r. The primary masses are given 
in lTetzlaff et al.l ()2010[ ). while the radii and primary star 
temperatures were estimated from spectral type rela- 
tions given in iCarroU fc Ostlid (|2006f ) . An eccentric orbit 
could have a larger maximum velocity than that esti- 
mated by equation [U if the orbit was oriented such that 
the the secondary star was moving directly towards or 
away from earth at or near periastron. A star in an ec- 
centric orbit cannot get so close to the primary that they 
touch, and its average distance must still be far enough 
to allow for the observed secondary star temperature. 
These conditions lead to a maximum eccentricity, given 

by 



1 



^inax(-^pi'im Rscc) 
G'(MpHm + Msec) 



(2) 



where i?prim and i?si 



are the radii of the primary and 
secondary stars, respectively, and w„iax is the maximum 
circular velocity given by equation [TJ Typical values give 
Cmax ~ 0.6. Eccentricities near this value could allow for 
velocities significantly greater than Umax given by equa- 
tion [T] 

The above analysis is summarized in Table |3l There 
are two binary candidate systems that we have not been 



5.1. HIP 26713 

The cross-correlation for this candidate is shown in fig- 
ure 01 The strong peak at -220 km s^^ has a significance 
just over 4cr, and corresponds to a 5600 K star model. A 
sensitivity analysis (Table [31 step 2) gives a median peak 
significance of ^ 6a with a large (^ 2a) spread. Due to 
the large spread, we cannot reject the peak based on the 
observed significance. 

The candidate radial velocity amplitude of 220 km s~^ 
is very near the upper limit given by equation [l] of 256 
km s"-^. If this candidate is a real binary companion, 
it must have been observed very near its radial veloc- 
ity maximum and the orbital inclination must be very 
near edge-on. Assuming a circular orbit and equal prob- 
abilities of observing any given phase or inclination, the 
probability of this occuring for this system is ~ 0.01. 
However, the probability for one system in our entire 
sample to be caught with this chance alignment rises to 
0.14. In addition, we cannot discount the possibility of 
an eccentric orbit leading to a higher value of Umax than 
estimated by equation[T] Without more data, we can nei- 
ther confirm nor disprove the existence of a companion 
star orbiting HIP 26713. 

If HIP 26713 is a true bina ry system, evolutionary 
tracks by iLandin et al.l (|2008[ ) give a secondary star 
mass of 1.6 ± 0.2 Mp,. The mass of the primary star is 
9.4±O.2M0 (Tetzl affet aIl l2010D. giving a mass-ratio of 
q = 0.17 ± 0.02. For a circular orbit, the corresponds 
to an orbital period of 10.0 days. If the companion star 
is on an eccentric orbit, its true period would be longer 
than this. 

5.2. HIP 92855 

There were seven observations of HIP 92855 on dif- 
ferent dates, all with the 2336 nm wavelength setting. 
The cross-correlations for four of the observation dates 
show a single strong peak when using a 6100K model 
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Fig. 5. — Cross-correlations for HIP 92855, for all dates observed. A 6100 K star model spectrum is used as the template for each 
cross-correlation. The y-axis is in units of the standard deviation of the cross-correlation function. A single strong peak is seen in the 
cross-correlations from 2007 May 9, 2007 August 2, 2008 September 19, and 2008 October 10. The reasonably strong peak on 2007 June 9 
is identified as probably arising from imperfect telluric line removal or random noise, since it has v > I'max given by equation [T] 



star as template. Figure [5] shows the cross-correlations 
of the telluric-corrected spectra with a 6100K model for 
all of the observations. Table |4] summarizes the sensitiv- 
ity and cross-correlation significance. The detection rate 
is the fraction of the 17 radial velocities between -400 
and 400 km that were correctly detected in the sen- 
sitivity analysis (Tabled step 2). The variation in detec- 
tion rate is due to the different signal-to-noise levels and 
telluric line corrections in the observations at different 
dates. The expected significance is the median signifi- 
cance of the radial velocities which were detected in the 
sensitivity analysis, in units of the standard deviation of 
the cross-correlation function. The observed significance 



and velocity are for the observed peaks. The velocities 
in Table |4] are corrected for the barycentric motion and 
the known systematic radial velocity of HIP 92855, while 
those in Figure [S] are not. 

As Table H shows, a 6100K star orbiting HIP 92855 
is at the limit of detectability with the direct spectral 
detection method. With the exception of the observa- 
tion on 2007 September 16, the cross-correlations for the 
dates with the highest detection rates have a single large 
peak. We consider this an excellent candidate for follow- 
up observations. 

If HIP 92855 is a true bina ry system, evolutionary 
tracks by iLandin et al.l (|2008l ) give a secondary star 
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Fig. 6. — Completeness diagram for the full sample of main-sequence B stars, split up by the effective temperature of the secondary star. 
The points correspond to the individual stars in the sample, and their sizes reflect the signal-to-noise ratio in the spectrum. Note that the 
signal-to-noise is calculated after the telluric line removal, and counts any telluric residuals as noise. The figures are also color-coded by 
the fraction of trials that detected the companion (see section [G]!. Contours are drawn to guide the eye. The red areas in each plot indicate 
the regions for wrhich our sample is complete. 



mass of 1.2 ± 0.2M(t,. The mass of the primary star is 
7.8 ± 0.2Mq (|Tetzlaff et al.ll2010n . giving a mass ratio of 
q = 0.15±0.04. The maximum radial velocity, observed 
on 2007 August 2, was 234 km s^^. If this is the radial 
velocity semi-amplitude and if the companion is on a cir- 
cular orbit, the binary orbit would have a period of 6.8 
days. If the true velocity semi- amplitude is larger cither 
because no observation was taken when the companion 
star was at quadrature or because the orbit is inclined, 
the period would be shorter than this. If the companion 
star is on an eccentric orbit, the period could be longer 
than 6.8 days if the large velocity observed was near pe- 
riastron. 

6. COMPLETENESS 

We now estimate the completeness of the direct spec- 
tral detection method applied to this data set. For each 
telluric-corrected observation, we created a series of syn- 
thetic binary-star spectra by adding stellar models to the 
data at various flux ratios, temperatures, and radial ve- 
locitie s. We used evolutionary tracks from iLandin et all 
(|2008f ) to find the luminosity of model companion stars 
with temperatures ranging from 3000K to 7000K in steps 
of 500 K. To find the model flux ratio Fg/Fp, we used 
the best-fit age for each star quoted in Table [TJ as well 
as the best-fit a.ge±aage- Model secondary spectra gen- 
erated with the Phoenix code (jHauschildt et al] I1999D 
were added to the telluric-corrected observations at 17 



different radial velocities ranging from -400 to +400 km 
s^^. Changing the radial velocity of the model spectrum 
changes where the companion spectral lines fall with re- 
spect to the telluric lines. Finally, we cross-correlated 
each synthetic spectrum with its corresponding model 
secondary star spectrum, and examined the cross corre- 
lation function. 

If the highest peak in the cross-correlation function 
was at the correct velocity, the companion was consid- 
ered detected. We then tabulated how many times the 
companion star was detected in the 17 radial velocity tri- 
als. Figure[6]shows the fraction of trials that detected the 
companion for all of the model radial velocities, as a func- 
tion of primary (B-star) mass and the binary mass-ratio. 
The points correspond to individual spectra, with their 
sizes indicating the signal-to-noise ratio in the spectrum, 
and the contours are drawn by interpolating between the 
points. The four panels are for different companion star 
temperatures. Companion stars with effective tempera- 
tures from 4600 — 5400 K have large regions with a very 
high detection rate. Stars cooler than about 4600 K are 
too dim to detect without much higher signal-to-noise 
ratios than present in our dataset, and stars hotter than 
about 5400 K do not have a strong CO bandhead and so 
the cross-correlation function is not as sensitive. 

Figure [6] shows that the direct spectral detection 
method is able to find companion stars with a mass-ratio 
of g ss 0.1 — 0.2, for a range of effective temperatures. 
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The regions with a detection rate near 1 are completely 
sampled, and a companion star in that region would be 
detected. For primary stars with M < IOM0, we are sen- 
sitive to almost all companions with 4600 < T^g < 5400 
K. 

7. MULTIPLICITY FRACTION 

We have not found any unambiguous low mass-ratio 
companions in our sample, though we do have two can- 
didates that require follow-up observations (HIP 92855 
amd HIP 26713). From the work described in section 
[6l we define a range of mass-ratios for which the direct 
spectral detection method is sensitive for each primary 
(B-) star. We can then rule out any companions with 
mass-ratios in that range, for that primary star. 

In order to convert these star-by-star limits on the pres- 
ence of a companion into upper limits on the multiplic- 
ity fraction of the parent population, we first count the 
number of stars that rule out companions in a particular 
range of mass-ratios. We then apply binomial statistics, 
where the probability P of finding k companions from n 
samples of a parent population with a true binary frac- 
tion p is given by 

mp,n)^^^Ai~Pr-^ (3) 

For no detected binary companions (k=0), the corre- 
sponding likelihood function for the binary fraction is 

P(p|fc = 0,n) = (n + 1) • (1 -p)" (4) 
A 90% upper limit is given by 

0.9= r {n + l)-{l-pYdp' (5) 
^0 

with solution 

P90 = 1-0.1^ (6) 

A similar derivation gives 90% upper limits for one de- 
tection (k=l). Figure [7] shows the 90% upper limits to 
the binary fraction as a function of mass-ratio. To find 
n in equation [51 we counted the number of stars that 
ruled out a companion at the given mass-ratio. We only 
counted companions that were found in all 17 radial ve- 
locity trials (see section |6l), and were always found with 
at least Aa significance. We also included upper limits 
assuming that HIP 92855, the more likely of our two can- 
didates, is a true binary systen i. Figure [3 a, l so sho ws the 
lower multiplicity limit set by [Evans et alj (|2011[ ) (blue 
dotted line) and the intrinsic O st ar ma ss-ratio d is tribu - 
tion derived by ISana et all ()2012f ). lEvans et al.l ()2011[ ) 
do not split their multiplicity by mass-ratio, and so the 
line shown in figure [7] is an average value. While we 
include them for comparison with our results, neither 
of these studies a re directly comparable to our sample. 
ISana et all (j2012D have only O stars in their sample and 
are only able to measure mass-ratios g « 0.2 — 1, although 
they consider mass-ratios down to q = 0.1 wh en deriving 
the intrinsic distribution. lEvans et al.l ()2011l ) are sensi- 
tive to similar mass-ratios as our sample, but they sample 
late B-stars (B4-B9) while we sample early B stars (BO- 
B5). Our results are almost perfectly complementary to 



those of lEvans et al.l (|20lH ). It is encouraging that our 
upper limits for 0.1 < g < 0.2, where our sample is most 
complete, lie in between the results of ISana et al.l (|2012D 
who measu re the binary fract ion of more massive pri- 
maries, and lEvans et al.l (|2011[ ) who measure the binary 
fraction of less massive primaries than we do. 

The above analysis assumes that the sample of B-stars 
given in Table[T]is representative of the B-star population 
as a whole. Most of the sample stars are field B-stars, 
which have a lower overall multiplicity than cluster or as- 
sociation B-stars (Mason et al. 2009). However, the close 
binaries this method is sensitive to would be difficult to 
disrupt with dynamical interactions in a cluster environ- 
ment, and we therefore expect this sample to be repre- 
sentative of both populations. A potential complication 
is that telluric standard stars are often chosen specifi- 
cally because they are not known binary systems. While 
this may introduce a bias for large mass-ratios, the direct 
spectral detection method is sensitive to low mass-ratio 
companions that may not have been found with other 
methods such as classical spectroscopy or imaging and 
so companions with q < 0.2 may be less effected. The 
amount of bias introduced into the above measurement 
depends on how carefully the individual observers chose 
the telluric standard stars, and so is very difficult to as- 
sess. 

8. CONCLUSION 

We have described a new technique for finding binary 
systems with a flux ratio of Fp/Fg « 100, where Fp and 
Fs are the fluxes from the primary and secondary star, 
respectively. In this technique, which we call the direct 
spectral detection technique, we use high signal-to-noise, 
high resolution spectra of a binary candidate. We remove 
the contamination from the Earth's atmosphere with the 
telluric modeling code LBLRTM, and cross-correlate the 
residuals with a library of stellar models for late type 
stars (F2-M5). A binary detection would appear as a 
strong peak in the cross-correlation function. 

We prove the feasibility of the direct spectral detec- 
tion method by adding a synthetic signal to real data, 
and successfully recovering it. We further investigate the 
completeness of the method in section [HI This method is 
sensitive to detecting a range of companion stars, set by 
the spectral type and age of the primary and the signal- 
to-noise ratio of the observation. Our sample is sensitive 
to almost all companion stars with 4600 < T < 5400, 
corresponding to binary mass-ratios of 0.1 < g < 0.2. 

We have applied this technique to a sample of 34 
archived main sequence early B-stars (B0-B5) with spec- 
tra taken with the CRIRES near-infrared spectrograph. 
We found no unambiguous companions in our sample, 
but identify two targets as candidate binary systems: 
HIP 92855 and HIP 26713. HIP 92855 is B2.5V type 
star, with a candidate companion star with effective tem- 
perature T = 6100 K and mass 1.2 ± O.2M0. Such a 
companion star is very near the detection limit of the 
direct spectral detection technique and deserves further 
follow-up observations. HIP 26713 is a B1.5V type star, 
with a candidate companion star with T = 5600 K and 
M = 1.6 ± O.2M0. This star was only observed once in 
our sample, and so may be a series of incompletely re- 
moved telluric absorption lines masquerading as a com- 
panion star. 
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Fig. 7. — Estimates of the binary fraction of B0-B5 stars, as a function of binary mass-ratio. This work found no unambiguous companions, 
and so we give 90% upper limits (soUd black line). 90% upper limits are also given assuming that HIP 92855 is a real binary system (dotted 
blac k line ) . The upper limits are only different within the la error bars on the mass-ratio for HIP 92855. The ne arly flat distribution found 
bv .Sana et al. (2012) is show n as the dashed red line. The average binary fraction found bv lEvans et al.l 11201 1) is also shown (dash-dot 
blue line). The i Evans et al.l 11201 II ) value is a lower limit and an average over all mass-ratios from 0.1 < q < 0.3, but they estimate that 
their sample is very complete, and so the true multiplicity fraction is quite close to their value. 

not expect a strong peak in the observed mass-ratio 
distribution. This interpretation may be supported 
by the flatness of the mass-ratio distribution as weU 
as disk instability si mulations that end with mas- 
sive c ompanions (e.g. IKrumholz et^\2m^ IClarkd 
[20091 ). 



We set upper limits on the binary fraction of early B 
stars as a function of binary mass-ratio (see Figure[7|) . As 
well as showing the upper limit for no detections in our 
sample, we also show the binary fraction upper limit as- 
suming the HIP 92855 is a true binary system. Our upper 
limits are strongest for mass-ratios q ~ 0.1 — 0.15, and are 
about 20%. We compare our limits to the intrinsic binary 
mass-ratio distrib ution for 0-type primaries derived by 
iSana et all (|2012f ). as well as the l ower l imit average bi- 
nary fraction seen bv lEvans et al.l (|2011[) for late B-stars 
(B4-B9). Our strongest upper limits (0.1 < q < 0.15) 
fall in between these two previous studies. 

Companion stars formed by circumstellar disk insta- 
bilit y are expected to have typical mass-ratios near q = 
0.1 feratter fc Matzne'3l2006l: IStamatellos fc WhitworthI 
[2oTiir near where our upper limits are strongest. If there 
was a large population of low mass-ratio companions 
formed by disk fragmentation, we would expect to see 
a peak in the mass-ratio distribution near g « 0.1. Since 
our results agree we U with the nearly flat mass-ratio dis- 
tribution derived bv lSana et al.l (|2012l ). it is unlikely that 
such a peak exists. There are several possible interpre- 
tations of this result, three of which we give below. 

1. Stellar companions formed by disk instability have 
a much lower characteristic mass-ratio than g = 0.1 
which remain invisible to observations. 

2. The mass-ratio distribution of companions formed 
by disk instability is very broad, and so we would 



3. Disk instability is not a dominant formation mecha- 
nism for low mass-ratio binary systems, and molec- 
ular core fragmentation alone can generate the 
nearly flat mass-ratio distribution down to low 
mass-ratios. 

It is difficult to distinguish between these interpreta- 
tions at this time. More observational work, as well as 
computational work, is required to explain the binary 
properties of high-mass stars. 

We would like to acknowledge Andreas Seifahrt, Rob 
Robinson, and Daniel Jaffe for their generous help with 
the telluric modeling procedure and some of the statisti- 
cal aspects of this work. We would also like to thank the 
referee for a quick review and many helpful comments. 
This research has made use of the following online re- 
sources: the SIMBAD database and the VizieR catalogue 
access tool at CDS, Strasbourg, France, and the ESO Sci- 
ence Archive Facility. Funding for this work was provided 
by a National Science Foundation CAREER award to 
Sarah Dodson-Robinson (AST-1055910) and by start-up 
funding from the University of Texas College of Natural 
Sciences. 



TABLE 1 Full star sample 



Star 


Spectral Type 


Age (Myr) 


Distance (pc) 


Maximum 
Separation (AU) 


Number of 
Observations 


HIP 23364 


B3V 


31.6 ±0.6 


31.65 


25.32 


1 


HIP 26713 


B1.5V 


7.2 ±2.5 


138.89 


111.11 


1 


HIP 27204 


BlIV/V 


12.6 ±4.6 


408.2 


326.5 


1 


HIP 30122 


B2.5V 


32 ±0.4 


111.1 


88.9 


4 


HIP 32292 


B2V 


8.2 ±0.1 


1111.1 


888.9 


1 


HIP 39866 


B3V 


25.1 ±2.6 


840.3 


672.3 


1 


HIP 48782 


B3V 


32.3 ±0.6 


370.4 


296.3 


1 


HIP 52370 


B3V 


17.2 ± 1.3 


58.14 


46.51 


2 


HIP 52419 


BOVp 


4 ±0.7 


250 


200 


2 
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TABLE 1 - Continued 
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30.12 


24.1 




HTP 7fil 9fi 
nil 1 uJ-^u 


B3V 


15.9 di 1.3 


62.89 


50.31 


]^ 


HIP 78820 


B0.5V 


13.8 lb 0.4 


123.9 


99.1 


1 


HIP 80582 


B4V 


50.1 ± 14 


19.96 


15.97 


2 


HIP 80815 


B3V 


10.5 lb 2.1 


95.24 


76.19 


4 


HIP 81266 


B0.2V 


5.7 lb 1 


175.44 


140.35 


12 


HIP 82514 


B1.5Vp+ 


20 lb 2 


50 


40 


1 


HIP 87314 


B2/B3Vnn 


23.2 lb 2.9 


43.1 


34.48 


7 


HIP 92855 


B2.5V 


31.4 lb 0.4 


31.85 


25.48 


7 


HIP 92989 


B3V 


7.9 lb 2.1 


126.58 


101.27 


1 


HIP 97611 


B5V 


45 lb 10 


66.67 


53.33 


1 



TABLE 2 Summary of adjusted line strengths. The units of line strength are 
cm~^/ (molecule X cm~^) 



Wavelength (nm] 


Molecule 


Old Strength 


New strength 


2317.12 


CH4 


5.445x10- 


-ai 


5.034x10- 


-•21 


2318.24 


H2O 


1.400x10" 


-24 


2.256x10" 


-24 


2328.51 


CH4 


2.521x10" 


-21 


2.371x10" 


-21 


2328.56 


CH4 


1.270x10" 


-21 


1.358x10" 


-21 


2340.12 


CH4 


3.085x10" 


-21 


2.963x10" 


-21 


2340.36 


CH4 


3.343x10" 


-21 


3.211x10" 


-21 


2351.64 


H2O 


1.670x10" 


-23 


1.393x10" 


-23 


2351.69 


H2O 


1.085x10" 


-23 


7.985x10" 


-24 


2352.43 


CH4 


3.144x10" 


-23 


4.031x10" 


-24 


2352.45 


H2O 


4.639x10" 


-23 


4.939x10" 


-23 


2353.62 


CH4 


2.708x10" 


-21 


2.654x10" 


-21 


2355.82 


CH4 


5.101x10" 


-21 


4.949x10" 


-21 


2358.9 


CH4 


5.160x10" 


-21 


4.710x10" 


-21 


2364.03 


H2O 


1.408x10" 


-23 


1.217x10" 


-23 


2367.23 


H2O 


2.078x10" 


-23 


2.182x10" 


-23 


2370.35 


CH4 


4.028x10" 


-21 


3.625x10" 


-21 


2370.41 


CH4 


2.437x10" 


-21 


2.021x10" 


-21 


2370.75 


CH4 


1.466x10" 


-21 


9.138x10" 


-22 


2371.39 


H2O 


3.905x10" 


-23 


3.171x10" 


-23 


236.62 


H2O 


1.146x10" 


-23 


9.186x10" 


-24 


2376.63 


H2O 


3.824x10" 


-24 


3.820x10" 


-24 


2378.2 


H2O 


1.134x10" 


-22 


1.021x10" 


-22 


2379.67 


H2O 


6.334x10" 


-24 


8.408x10" 


-24 


2385.98 


H2O 


6.051x10" 


-23 


5.407x10"^ 


23 



TABLE 3 Summary of Cross-Correlation Function Peak Rejection Steps. 
See Section [5] for more information 



Step 


Description 


Method 


1 


Telluric Residual Peak Identification 


Compare cross-correlation function of telluric- 
corrected spectrum with that of an uncor- 
rected spectrum. 


2 


Sensitivity Analysis 


Check that signal-to-noise ratio is high enough 
to detect secondary candidate at a range of 
velocities 


3 


Velocity Analysis 


Check that a blackbody with the candidate 
temperature can exist as close to the primary 
B-star as the velocity indicates (assumes cir- 
cular orbit) 



Detection of Low Mass-Ratio Stellar Binary Systems 



13 



TABLE 4 Summary of HIP 92855 observations. Significance is in units of 
the standard deviation of the cross-correlation function. The detection rate 
and expected significances are from the sensitivity analysis (Table [S] step 2). 



Date 


Detection Rate 


Expected Significance 


Observed Significance 


Velocity (km s '-) 


2007 May 9 


0.35 


3.3 cr 


4.1 cr 


-74 


2007 June 9 


0.18 


3.6 a 


N/A 


N/A 


2007 July 25 


0.35 


3.7 a 


N/A 


N/A 


2007 August 2 


0.47 


3.6 cr 


3.5 cr 


-234 


2007 September 16 


0.82 


4.2 a 


N/A 


N/A 


2008 September 19 


0.71 


3.4 a 


4.1 cr 


-126 


2008 September 10 


0.59 


3.3 a 


4.0 cr 


116 
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